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Abstract 

We study a rapidly rotating Bose-Einstein condensate confined to a finite trap in the framework of 
two-dimensional Gross-Pitaevskii theory in the strong coupling (Thomas- Fermi) limit. Denoting the 
coupling parameter by f/e 2 and the rotational velocity by Q, we evaluate exactly the next to leading 
order contribution to the ground state energy in the parameter regime loge| < n < l/(e 2 |loge|) 
with e — > 0. While the TF energy includes only the contribution of the centrifugal forces the next 
order corresponds to a lattice of vortices whose density is proportional to the rotational velocity. 

MSC: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q. 
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1 Introduction 

Bose-Einstein condensates respond to rotational motion of the enclosing container by the creation of 
quantized vortices. This remarkable manifestation of superfluidity has been studied, both experimentally 
and theoretically, in dilute, ultracold Bose gases since almost a decade and still offers a number of 
unsolved problems. We refer to the monograph [T], the review article [2J, as well as the papers [5]-[T7] 
for extensive lists of references. Most of the theoretical work has been carried out within the framework 
of the Gross-Pitaevskii (GP) equation for the wave function of the condensate. In the GP equation 
the interaction is encoded in a single parameter g = AirNa/L, where a is the scattering length of the 
interaction potential, N the particle number and L the length scale associated with the external confining 
potential. For rotating gases in their ground state the GP equation was derived in |10] (upper bound) and 
(lower bound) from the quantum mechanical many-body Hamiltonian with purely repulsive, short 
range interactions and fixed values of the rotational velocity and the coupling parameter as N — ► oo. 
The extension of this derivation to the case when the coupling parameter and the rotational velocity 
tend to infinity (or approach a critical value in the case of harmonic traps) has not yet been completed, 
but the leading order asymptotics of the many-body energy for large coupling and rotational velocity in 
anharmonic traps was computed in 12]. 
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Detailed results on the emergence of vortices as the rotational velocity is increased have been obtained 
within two-dimensional GP theory when the GP interaction parameter is large ('Thomas-Fermi' limit) 
and the rotational velocity is of the order of the logarithm of this parameter [HI HI IS EB] • In this case 
the number of vortices remains finite as the interaction parameter tends to infinity. The rotation has no 
effect on the energy to leading order in the coupling parameter but there is a logarithmic contribution 
due to the vortices in the next to leading order. By contrast, the papers |13j-[17j are mainly concerned 
with the situation when the rotation is so fast that the centrifugal energy and the interaction energy are 
comparable in magnitude. This holds when the rotational velocity increases like the square root of the 
interaction parameter. This case was also considered in the many-body context in |12j . relying partly on 
estimates from [16 . The main result of [16] and [17] was the rigorous evaluation of the GP ground state 
energy to leading order in the interaction parameter in the regime just mentioned. This energy can be 
computed by minimizing a simple density functional that contains besides the interaction term another 
term of the same order corresponding to the centrifugal potential. This functional was first introduced 
in |13] where many of the basic insights about the physics of rapidly rotating condensates in anharmonic 
traps can be found, see also [H] and [15] . We note that, since the rotational velocity is unbounded, the 
confining potential must increase more rapidly than quadratically with the distance form the rotational 
axis in order that the centrifugal forces do not tear the condensate apart. 

In the present paper we evaluate exactly the next term in the asymptotic expansion beyond the leading 
contribution in the parameter regime 

|loge| I/(e 2 |log £ |) (1.1) 

where the coupling parameter g has been written as 1/e 2 with e — > 0, and is the rotational velocity. We 
remark that the dimcnsionless parameter e can be interpreted as the ratio between the 'healing length' 
(47t7V/L 3 ) -1 / 2 and the extension L of the confining trap. The subleading term in the energy corresponds 
to the energy of a lattice of vortices of degree one such that the total vorticity is proportional to the 
rotational velocity. In order to bring out the salient points as simply as possible we restrict ourselves to 
the model considered in [IB], i.e., the case of a flat, circular trap ('bucket') of finite radius. 

When computing the upper bound on the energy we make a variational ansatz with a wave function 
that is essentially the product of a shape function, taking the deformation due to the centrifugal forces 
into account, and a function corresponding to a lattice of vortices uniformly distributed over the trap. 
The evaluation of the energy can be cast in the form of an electrostatic problem with the vortices 
playing the role of point charges while the vector potential due to the rotation can be regarded as an 
electric field generated by a uniform charge distribution. The optimal arrangement of the vortices is then 
determined by a minimization problem for the total electrostatic energy. When the rotational velocity 
reaches 0(1/ (e 2 \ loge|)) a different trial function, with the vorticity concentrated in a region where the 
density is small ('giant vortex'), gives a lower energy. This transition was first noted in [13j and the 
estimates of the present paper corroborate it since our rigorous upper bound to the energy is smaller 
than the energy of the giant vortex if <§C l/(e 2 | loge|)). 

To prove the lower bound the problem is reformulated in such a way that results from Ginzburg-Landau 
(GL) theory obtained in |20j and [21] can be employed. The strong inhomogeneity of the density in fast 
rotating condensates causes problems that make the reduction to the GL case not entirely straightforward, 
but once these have been overcome a lower bound that matches the upper bound to subleading order in 
the asymptotic parameter range (|1.1[) can be derived. The techniques of [5D] and [H] also turn out to be 
useful for the investigation of the vorticity of the minimizcr. 

2 The Mathematical Setting 

We now recall the setting of [TBJ that will be used in the present paper. The condensate is confined to 
the two-dimensional unit disc B\ and the rotational axis is perpendicular to the disc and passes through 
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its center. We note that this model can also be applied to the description of a three-dimensional rotating 
condensate confined to a long cylinder. The plane of the disc is the xy plane and r = (x,y) is the position 
vector with length r, while e 2 denotes the unit vector in the z direction. The complex valued order 
parameter (the wave function of the condensate) is denoted by W(f). In the non-inertial rotating frame 
the GP energy functional can be written as 



rGP 



= J df I (V - iAj * 



(2.1) 



where the vector potential A is given by 

A=^e z /\r. (2.2) 

For convenience we also introduce the abbreviation 

Lj = eVL. (2.3) 

For fixed ui the centrifugal and the interaction terms in ()2. 1(1 are both 0(l/e 2 ). The kinetic first term, 
containing A ~ fi, is formally also of order 1/e 2 if il ~ 1/e, but a complex phase factor in due to 
vortices, can partly compensate the effect of A. Indeed, in the ground state this term is of lower order 
as we shall see. 

The ground state properties of the condensate are obtained by minimizing the GP functional over the 
domain 

D GP = {*G J ff 1 (S 1 )|||*|| 2 = l}. (2.4) 

Here H 1 (Bi) denotes the Sobolev space of complex valued functions * on B\ such that both * and V* 
are square integrable. The choice (|2.4|) naturally leads to (magnetic) Neumann boundary conditions for 
the minimizer on dB\. Alternatively one could impose Dirichlet boundary conditions. For Q ~ 1/e this 
would affect the energy to order 0(l/e) that is negligible compared to the vortex contribution <3(f2| loge|) 
that we are interested in. In other parameter regions the effect of the boundary conditions can be more 
significant, and the same remark applies to an extension of our analysis to homogeneous potentials as in 
|17j . For simplicity we shall, however, in this paper stick to the choice (|2.4j) that highlights the vortex 
contributions. 

We denote by E the GP ground state energy and by * GP any corresponding minimizer. The 
existence of such minimizer(s) as well as the fact that any minimizer solves the GP differential equation 



2 



-iA) V GP -A 2 V GP + 2s- 2 \y GP \ 2 V GP = f i GP V GP , (2.5) 

with boundary condition V r 5 ,GP = on dB\ can be deduced by standard techniques (see, e.g., (7|). The 
chemical potential p GP is fixed by the L 2 — normalization of * GP , i.e., 

/i GP = jB GP +e -2||^GP||4 (26) 

In [IB] we studied the asymptotics of E GP as e — > and proved that the energy is well approximated 
to leading order by minimizing the 'Thomas-Fermi' (TF) functional 



Note that, unlike in [T51 [T7], we have included the factor 1/e 2 in the definition of £ TF . The density 
p(r) > is the probability density associated with a condensate wave function *, i.e., p — 1*1 2 . The TF 
ground state energy, 

E TF = min £ TF [p}=£ TF [p TF ], (2.8) 



IIpIIi=i,p>o 
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and the corresponding normalized density p can be explicitly calculated. The formulas and some 
properties of relevance for this paper are collected in the Appendix. We note in particular that the 
centrifugal forces may create a 'hole' in p TF , i.e., the density p TF vanishes on a disc centered at the origin 
if ijj > L>h = 4/v / 7T- 



3 The Main Results 

The main results proved in are, in a slightly different notation, contained in the following 

Theorem 3.1 (Leading order ground state energy asymptotics [16j) 

7/f2 ~ 1/e as e — > 0, then 

E GP = E TF + 0(E- 1 \loge\), (3.1) 

whereas if 1/e <C tt, 

E GP = E TF + 0(£ - 2) + (( eJ2 )2| log£ |)_ (3.2) 

In this paper we investigate the correction beyond the leading order TF term for the parameter range 
| logs] <fl< l/(e 2 | loge|). Our main result is as follows (the notation Q < 1/e means that f2 < C/e as 
e — * 0, with some C < oo): 

Theorem 3.2 (Improved ground state energy asymptotics) 

If | loge| < Q, < e _1 as e — > 0, then 



whereas, if e 1 < fl <C e 2 |loge| 1 , 



£;^ =jE ;^ + ^^n(l + (l)). (3.4) 

In [H] we have shown that as a consequence of the energy asymptotics |\I/ GP | 2 converges as e — > to 
p TF in L'-norm. Inside the hole, if present, it is exponentially small, i.e., bounded by exp(— const. /e 13 ) 
for a (3 > 0. See [IB], Propositions 2.4 and 2.5. 

The energy bounds also allow to prove a result about the uniform distribution of the vorticity of v]/ GP 
outside the hole, at least for < 1/e: 

Theorem 3.3 (Uniform distribution of vorticity) 

Let <3> GP be any GP minimizer and e > sufficiently small. If | loge| <C Vl < e _1 , there exists a finite 
family of disjoint balls {B l £ } C supp (p TF ) such that 

1. the radius of any ball is smaller than J7~ 1//2 . 

2. the sum of all the radii is much smaller than J7 1 / 2 , 

3. on dB\, |* GP | > Cllogte 2 ^- 1 with C > 

and, denoting by ri^ £ the center of each ball B\ and by d^ e the winding number of |'3> GP |~ 1 v]/ GP on dB\, 

^ J2 di, s 6 (r- n, s ) X TF (r) dr, (3.5) 

in the sense of measures, where X TF (^) stands for the characteristic function o/supp (p TF ). 
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For Q ^> 1/e the vorticity distribution is still an open question. In this regime ty GP is not uniformly 
bounded in e and is essentially supported in an annulus of very small width ~ lo~ 1 close to the boundary. 
As we shall see, a trial function with a uniform distribution of vortices still gives the right energy to 
subleading order for <C l/(e 2 | logeQ, but for larger Q, a trial function without any vortices in the 
support of /5 TF (a 'giant vortex') has lower energy. It can be expected that for the true minimizers the 
vortices are gradually expelled from the essential support of the density as SI approaches l/(e 2 |loge|) 
but there are so far no rigorous results on the details of this phenomenon. The numerical investigations 
of [22j . however, support this picture. 



4 Energy Upper Bound 

For an upper bound we test the functional (|2.1[) with a trial function of the form 



^{r)=c^)i(f)g(r), (4.1) 

where c is a normalization constant, p a suitable regularization of p TF , g is a phase factor, and £ a 
function that vanishes at the vortices, i.e., the singularities of g. To define the functions precisely we first 
introduce some notation. 

We denote by £ a finite, regular lattice (triangular, rectangular or hexagonal) of points ft 6 By. Each 
lattice point f* lies at the center of a lattice cell Q l and the lattice constant £ is chosen so that the area 
of Q 4 is 

\Q' l \ = ^- (4.2) 

Thus, 

I = (const. )fT 1/2 (4.3) 
and the total number of lattice points in the unit disc is 

7V = |(i + o(n- 1 / 2 )). (4.4) 

For large u> — eil the support of p TF has an area of the order (uj + l) -1 and the number of lattice points 
on the support of p TF is of the order 

Af 1 = (w + 1)^. (4.5) 

In particular, for O > l/e, M 1 = 0(l/e). 

Using complex notation £ = x + iy for the points r = (x,y) GM 2 the phase factor g is defined as 

The phase factor is singular at the lattice points but these singularities are compensated by the function 

41 J I t^K-Cil if |C-Ci|<*. [ 

Here t, with 

min{e, (e/ft) 1 / 2 } < t <C fT 1 / 2 , (4.8) 

is a variational parameter that will be fixed later. Thus £(r) vanishes at the lattice points fi and is equal 
to 1 outside of the union of the discs B\ of radius t centered at those points. 

The size of t can be estimated by the following heuristic argument. The kinetic energy of a vortex in 
a cell is of the order ekin ~ //(l/ r ) 2r ^ r ~ ^>g(£/t). Creating a vortex also causes an excess interaction 
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energy because the density depletion in the vortex core of radius t has to be compensated by an increase 
in density elsewhere. This leads to the additional interaction energy e- ln t ~ plt/e) 2 , and by minimizing 
ekin + ei n t we obtain t ~ ep^ 1 ! 2 . For slow rotations where p = 0(1) this gives t ~ e, while for rapid 
rotation, where p — 0(u>), we obtain t ~ (e/il) 1 / 2 . These heuristic considerations are confirmed by the 
rigorous estimates below. 

The density p(f) can for ui < uj-h (see Eq. (jA.l|l ) simply be taken to be equal to the TF density 
p TF (f) (|A.3|) (note that p TF depends on oj). For ui > ut^, however, p TF vanishes in a 'hole' of radius 
i?h = 1 — const. (a;) -1 (see Eq. (|A.6[) ) and \J p TF does not have finite kinetic energy. Hence it is necessary 
in this case to regularize p TF near the boundary of the hole. At the same time one has to take care that 
the TF energy of the regularized density remains close to E TF . Both conditions are met if we put 

f if r<R h , 

p(r) = i p rF (R h + Q- 1 )n 2 (r-R il ) 2 if R h <r<R h + Q- 1 , (4.9) 
[ p TF (r) otherwise. 

Thus the regularized density is equal to p TF except in an annulus of thickness f2 _1 around the hole, where 
it increases quadratically with the distance from the hole. The latter property ensures finiteness of the 
kinetic energy. We also note that, by (|A.6|I and 1|A.7|) . p TF (-Rh + ^ _1 ) = 0(e 2 Q), so that, for any fe Bi, 

p(r) = p TF (r) + 0(e 2 n). (4.10) 

We now collect some simple estimates that are needed in the proof of the upper bound. In the 
following C will stand for a positive, finite constant that may differ from line to line but is independent 
of and e. 

First, note that p TF < C(uj + 1) while the area of the support of p TF is C(u> + The density of 

vortices is S1/27T and the area of each vortex disc is irt 2 . Also, e 2 il = o(l) by assumption. We thus have, 
using flUU) and P~TU]) . 

/ pf = / 9 - J P(l - C 2 ) > / P TF - e 2 - Cn ■ t 2 > 1 - 0(t 2 n). (4.11) 
Hence the normalization constant satisfies 

c<l + Ct 2 fl. (4.12) 

Likewise, using that |V£| = t^ 1 in each vortex disc and zero outside the union of the discs, while the 
number of vortices in the support of p is < (u> + l) -1 ^, 

\\y/P v ^ll2 ^ C ( UJ + 1) ' *~ 2 ' ( w + !r lfi -t 2 = Cn. (4.13) 
Next we consider, for w > Wh, i.e., Rh > 0, 

|Vr T1 ''- 

r<R t +n- 1 9 4 Jr^Rm+U- 1 

By (031) and ([ITUf first term is bounded by CO. ■ (e 2 fl). Using (jA3j) we obtain 



M2 ^[llf<lf M. + U 11^1. ,4,4, 

4 7 p 4y r<R +n -i p 4 J r>J i m+ [i-i p lt 



f \l^l = C ( £ nf f -^!^<c(^) 2 f ^^c^nj-niiogei. (4.15) 

T he ab ove estimate shows that the closer is to e~ 2 , the larger is the kinetic contribution of the profile 
yj p TF and for ft ~ (e 2 |loge|) _1 is becomes of the same order as the other remainders, i.e., ~ Q. For 
w < ^h, on the other hand, p = p TF , and V^/p is uniformly bounded in to, so H^V^/pHl < C in this case. 
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Because g is a phase factor while J~p and £ are real- valued functions we have 



(v-*a)(vp£s)| 2 = |v(Vp£)| 2 + £ 2 



V - iAj g < 

2£ 2 |VVp| 2 + 2 p |V£| 2 + £ 2 p 



v-u )g 



(4.16) 



We now obtain, using Eqs. (|4.12p ~ (|4.15p . 

£ GP [*] - £ TF m 2 } =c 2 I dr (V - il) (VP C 5) 



< 



dr£ 2 p 



V-iA U 



2c 2 



df £ 2 |VVp| 2 + 2c 2 / dr p |V<er < 

JBi 

2 

5 



2 



(l + 0(* 2 fi)) J dr£ 2 p TF (v-iA e ^g +C{n + e 2 n 2 \\oge\}. (4.17) 



The estimate for the vortex kinetic energy J B dr£ 2 p TF | (V — *^4) g\ 2 is given in the following Proposition 

Proposition 4.1 (Vortex kinetic energy) 

Ifs—yQ, and 1 < O < 1/e 2 , then 



2 „TF 



df£ 2 p 



(V - iXj 2 < Jfi | log(< 2 fi)| + 0(0.) + 0(fi (e 2 ft) 1/2 | log(t 2 ^)|). 



(4.18) 



Proof: The idea behind the proof is the electrostatic analogy that was already mentioned in the 
Introduction and is made precise in Eq. (|4.22p below and the considerations following it. The vortices, i.e., 
the singularities of the phase factor g, play the role of unit charges while the vector potential corresponds, 
after a conformal transformation, to the electric field of a uniform charge distribution. The density of 
the vortices is chosen in such a way that the field from the uniform charge distribution is compensated 
as far as possible and this requires in particular that each unit cell Q 1 has total charge zero. If the unit 
cells were rotationally symmetric there would be no interaction between them by Newton's theorem. 
Complete rotational symmetry is, of course, not possible, but the closest approximation to it among the 
regular lattices is a lattice with hexagonal cells, i.e., a triangular arrangement of the vortices, that gives 
the lowest electrostatic interaction energy. However, the difference between the three possible types of 
unit cells, triangular, rectangular and hexagonal does not show up in the term of order f2 1 log(i 2 S7)| but 
only in higher order corrections to this contribution. 

To formalize these ideas we note first that |(V — iA) g\ 2 = |V</> — A\ 2 where (f> = arg(£ — d) is the 
phase of g. The conjugate harmonic function 



E lo si 



satisfies 



where V r = e r ■ V = d/dr and V# = e$ • V = r~ 1 d/d'd. With A = A(r)e$ we thus have 



(4.19) 
(4.20) 



W<j>-A 



V, 



V r <j) - A 



e$ + V r e r — A e. t 



V0 - Ae r 



(4.21) 



We now define 



E(r) = V4>(r) - A(r)e r 



(4.22) 
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and note that V(j> = — fi) /\f — fi\ 2 can be regarded as the electric field generated by point charges 
localized at the positions of the vortices, while A(r)e r — (Sl/2)re r is the field generated by a uniform 
charge density of magnitude Sl/2ir = jQ 1 !" 1 . We can thus write E(r) — J^i Ei(r) — J2i V$j(r) with 



$i(r)= df a^F) log \r-f\ (4.23) 

and 

otf) = 6{? - n) - {QT'xi^) (4.24) 
where Xi is the characteristic function of the cell Q l . 

By a transformation of variables, writing r = SI -1 / 2 :?, we map the cells Q l of side length I ~ f2 -1 ' 2 onto 
cells Q\ of side length O(l). The characteristic function of Q\ is denoted by Xi,i (&) an d we use the index 
1 also for the charge densities, electric fields and potentials generated by the cells Q\. We can then write 

Ui{f) = Q [S^ - Xi) - IQir^i.iCf )] =fiff j) i(/) (4.25) 

and 

Bj(f) = fi 1/2 £!t,i(f) (4.26) 

where 

= V / dx' ^(f )log|z -x 1 1 = V*i,i(£). (4.27) 

JBi 

Consider now the cell Q\ centered at the origin. The multipole expansion of <&o,i {%) for x Ql is 

v-^ Cfc cos(fctf) + Sk sin(fc#) 
$o,i(*)=?log|^|-X) nk — ( 4 - 28 ) 



fc=i 



with 



9 = / df cr ,i(£')> 

Ck^k- 1 { dx' a Q i{x")\x'\ k cos(M'), Sk^k- 1 [ dx' a i(x')\x'\ k sm(M r ). (4.29) 

By neutrality of the charge distribution (|4.24p it is clear that q = and by symmetry of the unit cell 
it is also clear that there is no dipole moment, i.e., C\ = S± = 0. We conclude that $0,1 (x) decays at 
least as |a?| 2 and the corresponding field Eo t i(x) decays at least as \x\~ 3 . For square or hexagonal cells 
it decreases even faster. 

All cells Q % are obtained by translations and scaling from the cell Q\. From the considerations above 
(note, in particular, Eq. (|4.26[) ) we can thus conclude that if two of the original cells, Q l and Q 3 have 
distance 0{Sl~ 1 ^ 2 n) from each other, then the strength of the field Ej(f) for r € Q l is at most 0(Sl 1 ^ 2 n^ 3 ). 
Since, for a fixed cell Q l , there are at most 0(n) cells at distance 0{Sl~ 1 / 2 n) from it, we can estimate for 
re Q % 

E{r)-Ei{r)\ <^\E 3 {r)\ < const. Sl 1/2 ^ n ■ n~ 3 = 0(S1 1/2 ). (4.30) 

i¥* " 

Writing 

\E\ 2 = \E t \ 2 + 2{E - Ei) ■E i + \E- E t \ 2 (4.31) 
and using the simple bound Ei(f) < \r — , we conclude that for fe Q l 

\E(f)\ 2 < \E t {r)\ 2 + const.{Q 1/2 \r- r^ 1 + fi) (4.32) 
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,cn-' 

dr \Ei{r)\ 2 < const. / 



/ df\E(r)\ 2 - dr |£;(r)| 2 < const. / dr r (ft 1 / 2 ^ 1 + fi) = 0(1) (4.33) 



while 



dr £(f) 2 |£(r)| 2 - / df f (r) 2 |£i(r)| 2 < const. / dr r (r/i) 2 (f7 1/2 r^ 1 + 0) = 0((< 2 0) 1/2 ). (4.34) 



On the other hand, since Ei(r) < \r — fi\ 1 , 

,CI!' ' 

lQ*\Bl Jt 

and 



/ dr \Ei(r)\ 2 < 2ir / drr r~ 2 = tt| log(t 2 0)| + 0(1) (4.35) 



dr £,{r) 2 \Ei(r)\ 2 < 2ir / dr r (r(n/e) 1/2 ) 2 r- 2 = 0(1). (4.36) 
B\ 7 

Putting all the estimates above together we obtain 

f dfp TF (r)ar) 2 \E(r)\ 2 < (l + 0((t 2 n)^ 2 )) V sup p TF (rO (tt| log(i 2 0)| + 0(1)) . (4.37) 

JBi V ' • r6Q* 

It remains to estimate the Riemann approximation error 



(4.38) 



We use here that ||d / o TF /(ir|| 00 < C(ef2) 2 and that the number of cells Q % that intersect the support of 
p TF is bounded by Oe _1 (l + fT 1 / 2 ). Hence 

11 < CQ- 1 ■ n- 1 / 2 ^) 2 ■ e-\l + fl- 1 / 2 ) = C (e 2 0) 1 / 2 (l + fl- 1 / 2 ). (4.39) 

It now follows that the right hand side of (|4.37[) is bounded by 

(1 + Odt 2 ^ 1 ' 2 )) (1 + K) lOV 1 [A \og(t 2 tt)\ + 0(1)) = 

\Sl I log(£ 2 0)| + 0(0) + 0(n [e 2 Vl) 1/2 \ log(t 2 0)|). (4.40) 

□ 

To complete the proof of the upper bound we still need to estimate the difference between f TF [|\E'| 2 ] 
and E TF = £ TF {p TF ) and choose the radius t of the vortex discs. 
The TF functional is 



£ TF [|*| 2 ] =e- 2 J dr j|vl/| 4 - 



(eO) 2 r 2 |*| 2 



(4.41) 



We consider the two terms separately. For the nonlinear interaction term we use that c = 1 + 0(t 2 D,) 
and £ 2 p < p TF to obtain 



e- 2 



f dr\^ = e- 2 f df(c 2 e P ) 2 < 1 + C / n f dr(p TF ) 2 =e- 2 f dr (p TF ) 2 +rcmainder. (4.42) 

JB, JB, £ J B\ JB, 
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Since p TF < C (eCl + 1) and J p TF = 1, the remainder is 



^/ df( P TF ) 2 <c5Wl) (4-43) 



and this has to be small compared to 57| log(i 2 57)|. If 57 < 1/e this is clearly satisfied for t — e. On the 
other hand, if efi 3> 1 we can take t = (e/57) 1 / 2 Note that for 57 ~ 1/e both choices coincide. 

The centrifugal contribution can by partial integration and using the normalization of $ be written 

57 2 r tt57 2 r 1 
/ drr 2 |*| 2 = + tt57 2 / drr$(r), (4.44) 



with 



Likewise, 



with 



$(r) = ( dr'r' \^{r')\ 2 . (4.45) 
Jo 

^ j drr 2 p TF = + tt57 2 j dr r $ TF (r), (4.46) 



$ TF (r) - / drVp TF (r') 



From (|4.9j) and (|4.12[) we obtain 

$(r) < $ TF (r) + C*t 2 57. (4.47) 
Moreover, the support of $ as well as 4> TF has area < (e57 + l) -1 . Hence 



57 2 / dr |$ TF (r) - $(r)} < C57 2 ■ t 2 57(e57 + l)" 1 . 
•A) 



(4.48) 



If e57 is bounded and t = e, this is bounded by C 57 • (e57) 2 < C57. If ef2 3> 1, we take i 2 = e/57 and 

obtain again C 51 as bound. 

We summarize the findings in the following 

Proposition 4.2 (Energy upper bound) 

For e — > and 1 <C 57 < 1/e we have 

e gp < ^tf + i Q | i og ( £ 2 )| + (n), (4.49) 

and for 1/e < Q < 1/e 2 

£ GP < £ TF + Ifi| loge| + 0(fl) + 0(57 (e 2 57) 1/2 | loge|). (4.50) 

As we will see in the next section, the upper bounds are matched by corresponding lower bounds 
only in the parameter range | loge| <C 57 <C l/(e 2 | loge|). In fact, for 51 < | loge| there are only finitely 
many vortices @ HH HE] and the upper bound (|4.49p is too large. For l/(e 2 |loge|) < 57 -C 1/e 2 , on 
the other hand, a trial function different from (|4.1| gives lower energy than (I4.50| . This is the trial 
function considered in [16] Eq. (3.36) for 57 3> e _1 that corresponds to a 'giant vortex' where all the 
vorticity is concentrated at the center and the support of p TF is vortex free. In fact, for such a trial 
function the next correction to the TF energy is 0(l/e 2 ), cf. (|3.2[) . and this is smaller than 57 1 loge| in the 
parameter range. This transition at 57 ~ l/(e 2 | logeQ can also be understood by the following heuristic 
argument, employing the electrostatic analogy: For 57 3> 1/e the number of cells in the support of p TF is 
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~ 1/e. Without vortices each cell has unit 'charge', originating from the vector potential, and the mutual 
interaction energy of the cells is of the order 1/e 2 . Putting a vortex in each cell neutralizes the charge so 
that the interaction energy becomes negligible, but instead there is an energy cost of order Q\ loge| due 
to the vortices. Equating these two energies leads to ~ l/(e 2 | loge|) as the limiting rotational velocity 
above which the ansatz (|4.1|) is definitely not optimal. 

It should be noted that also for 1/e < Q <C l/(e 2 | log e | ) one could for the upper bound replace the 
distribution of the vorticity on a lattice within the 'hole' by a single phase factor corresponding to a giant 
vortex at the origin, but in order to obtain the correction beyond the TF term the support of p TF can 
not be vortex free. The detailed vortex distribution of the true minimizer of the GP energy functional is, 
however, an open question. 



5 Energy Lower Bound 

The lower bound to the GP ground state energy E GF will be proved by a step-by-step reduction to the 
lower bound of the energy of a Ginzburg-Landau (GL) energy function for which results of [HH [2H1 [H] 
can be employed. As a preparation we first prove a bound on the GP minimizers in terms of the TF 
density: 

Lemma 5.1 (Upper Bound for |\1/ GP |) 
For e -> and | loge| < f2 < (e 2 | loge|) _1 , 



|* GP f <P TF d)(l +o(D). 



(5.1) 



Proof: Setting U = |\1/ GP | , we first note that the upper bound (|4.50|) and the trivial lower bound 
E GF > E TF imply the convergence of U to p TF in L 2 -norm. Indeed, using the simple bound 2p TF > 

£ TF + £- 2 ||p TF || 2 , we have 



e 2 /i TF + u> 2 r 2 /4, the L 1 normalization of U, and the identity /x TF 



dr (U- p 



< 



dr 



U 2 -p TF U- 



Lu 2 r 2 U 



„TP' 



= e 2 (£ [U] - E 



<e 2 (E GF -E TF ) < Cne 2 \\oge\ = o(l), 
by (|4.49p and (|4. 50(1 and the conditions on fi. As a consequence 



\U\\ 



I TFII 2 o 
\P 2 = 2 



df P TF (u- P TF )+ / df ([/-p TF )^< (0 TF (l) 1 / 2 (l) 



(5.2) 
(5.3) 

(5.4) 



TF||2 



< 



where wc have used the Schwarz inequality and the trivial bound \\p 
follows from the L 1 normalization of p TF . Now p TF (l) > C(u> + 1) and theqrefore 



\P 



TF I 



p ijH (l), which 



Since e 2 (/i 



GP 



TF I 



I2 < (1)/"(1). 

o TF 1 1 we thus have 

^-D<o(l)Al). 
Now acting as in the proof of Proposition 2.4 in [TE\, we obtain from the variational equation 



re - 1 

TF ) = e 2 (E GF -E TF ) + \\U\\ 2 - ] \p ,,, 

2c. .GP ..TF^ ^ \ „TF 



(5.5) 
(5.6) 



\au< 



e 2 fi GF 



T~ 2U 



|<[e 2 (^ F - M TF ) + 2(p-(l)- C /)]|< 



U 



by (|A.5p and the above estimate for e 2 (p, GF 
nonnegative and thus (|5.ip holds. 



.TF") 



2[(l + o(l)) P TF (l)-C/]-, (5.7) 
At the maximum of U the left hand side of (|5.7[) is 
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□ 



We now proceed with the proof of the lower bound. The first step is the extraction of the TF profile 



p TF from the GP minimizer 5 ,GP , i.e., the ansatz *S> 



GP 



'u, which, one the one hand, allows to get 



rid of the leading order term in the energy asymptotics and, on the other hand, implies that u minimizes 
a weighted GL functional. Unfortunately such a factorization is well defined only if the TF profile p TF 
does not vanish inside B± , i.e., for u> < lo^. In order to get rid of this problem we first restrict the 
integration domain in the GP functional and set 



with 



by (fTTTj) . Note that 



T= {reBi \p TF (r)> cj\ log S]- 1 } 

S = e 2 n\\oge\ < 1, 
|log5i < Cjlogel, 



(5.8) 

(5.9) 
(5.10) 



since 5 > e 2 |loge| 2 , by dHU), so that > \ogS > log(e 2 | loge| 2 ) > Cloge. Note that Rl+uj- 1 ] log^" 1 < 
1, since R 2 = 1 — Cuo^ 1 and | log <5 1 ^> 1, so, by (|A.7|) . the set Tis not empty. Moreover, if ui/uiyy is 
sufficiently small then the set T coincides with the whole trap B\ since, in that case, p TF (r) > C > for 
any r € B\. 

For fl and e satisfying (II. ip we now define for r E T 



u(r) = y GP (r)p TF (r) 



-1/2 



(5.11) 



This is a smooth function with |u| < C| log (5| because of Lemma [5.11 Adding the kinetic energy term 
to both sides of (|5.2p we obtain, exploiting the nonnegativity of the integrand, 

E GP >E TF +f dr ( (V-iA)V GP 2 

Introducing the weighted GL-type functional 

2 



• e 



p TF I^GPI 



(5.12) 



£ GP [ 



u 



dr p 



TF 



T 



(r){|( 



iA) u 



e- 2 p TF 



(r)( 



1 



(5.13) 



we thus obtain, since Tc supp (p TF ), 



1 

2 .,, 



g GP ^GP] _ £TF _ £ ~gp [ u ]>t I cfrVp TF ■ V \u\ 2 > Coj 2 I drr- V \u\ 2 > 



Cuj 2 J df \u | 2 > -Cw|log S\, (5.14) 



which yields, by (|5.10|) . 

E GP >E TF +£ GP [u}~Cu\\oge\. (5.15) 

According to (j5 . 1 5[) the correction to the leading term E TF can thus be estimated from below by a 
weighted GL energy £ GP [u] , where the Lebesgue measure is replaced by p TF (r) df. Compared with the 
usual GL setting there are two differences, however: The internal magnetic field A is in our case fixed 
from the outset and the coupling parameter is p TF (r)e~ 2 , i.e., it depends on the TF density at each 
position. 
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To deal with the latter point we decompose the integration domain into small cells: Let C be the 
square regular lattice 

t = jfj = (ml, n£),m, neZ\Q l C t| , (5.16) 



where Q 1 denotes the lattice cell centered at f\ <E L and 



loge| 



n 



<C £ -C min 



1, 



1 



w\log6\ 



(5.17) 



Note that the above conditions are compatible: Multiplying both side by ui and assuming that w > 
| log 5 1 ~ 1 , (|5.17p becomes \/~5 <§; loI -C | log(5| _1 , which can always be fulfilled since (5 <C 1 by definition. 
Note also that the lattice spacing is much larger than the one chosen in the upper bound proof, where 
the lattice constant was I ~ SI -1 / 2 . By the lower bound on £ each lattice cell can be expected to contain 
a large number of vortices which turns out to be helpful for estimating the energy. The upper bound on 
I guarantees that I is much smaller than the width of T, which is of order (u> + 1) _1 . This is useful for 
the extraction of the TF profile. 

By (|A.3|) . p TF (r) > p TF (ri)(l — O(£oj\ log <5 1 ) ) , for any re Q 1 , so that the above inequalities imply 



£ GP [«] > E / . dr> TF (r) ( I (V - i A) 



+ £ -V F (r)(l-N 



> 



(i- (i))Ep TF (n)^)H, (5.18) 



with 



dr 



(v-iA)u\ 2 + e- 2 p TF ( n )(i-\u\ 2 y 



(5.19) 



Now the analogy with the GL functional is made explicit, since, except for the coupling parameter which 
still contains p TF , the functional £^ 1 ' is precisely the GL energy functional 



£ GL [u,A'] = / df 



V-iA' 



+ 



VAi'-L 



V TF (r.) (l-H 2 ) ; 



evaluated at (it, A) with an external magnetic field h cx — Qe z . It is clear that the GL energy 

E GL = inf £ GL [u,A'} 

u,A' 



(5.20) 



(5.21) 



is a lower bound to the ground state energy of £ ^' because the configuration with the uniform internal 
magnetic field corresponding to A' = A = ile z A r/2 is only one among all possible configurations 
considered in the minimization of the GL functional. 

We can now state the main estimate needed for the proof of the lower bound. 

Proposition 5.1 (Lower bound inside cells) 

For any satisfying (|1.1[) and e sufficiently small, it is possible to find £ in such a way that (|5.17[) is 
fulfilled and 

£M [„] > f*!£l2l (1- (1)), (5.22) 



where 7 = min[e, e 2 tt] . 
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Proof: The key point in the proof of (|5.22[) is a rescaling of Q l (together with the choice (|5 . 1 7[) of the 
lattice spacing), which allows to reduce the problem to the minimization of a GL functional in a different 
regime. We thus set x = £ (f — fj), 



u(x) = u (?i + ixj , B(x) = I A (?i + ix\ , 
where f% stands for the center of Q l . By such a change of coordinates in (|5.19p , we obtain 

*» H - *•> w - / {| (v - ,s) «|" + (i - ,,f } , 



(5.23) 



(5.24) 



where Q.% is a unitary square centered at the origin. Note that the rescaled vector potential B is explicitly 
given by 



B(x) 

and the corresponding magnetic field is 



flle z A r, fl£ 2 e z A f 

7. ~t~ ~ i 



h = curlB = nl 2 . 



(5.25) 



(5.26) 



In the following we investigate the minimization of the functional £^>: We first notice that, by gauge 
invariance, one can get rid of the constant term fl£e z A f\/2 in (|5.25p : 



inf £ {i) [u] > _ inf 



da; 



We now introduce a new infinitesimal parameter e defined as 



(V - i£ 2 A{x)) uf + e- 2 £ 2 P TF (n) (l - \u\ 2 )' 



(5.27) 



by 



and (f5~T7]) . It follows that 



< C7e« 



' n|log<S| 
w| log el 



< C7e« 1, 



(5.28) 



^'[u]> inf / df 



V ■ 



i/i ex e, A £ 



l-M 2 



(5.29) 



for a magnetic field /i ox satisfying the conditions 



|ioge| < /i cx = ni 2 < — . 



(5.30) 



Indeed, by (SHJ> and (jOS]) . 



fir = 



e 2 P TF (rO 



<^«1, ^ 2 »|log £ |>|log £ | 



(5.31) 



because > loge = loge — log{£ \J p T¥ '(r^)) and 



/° TF ('"i) <C min 



«1, 



(5.32) 
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which implies loge > loge and | loge| < | loge|. 

The functional on the right hand side of (|5.29[) is precisely the GL functional on Qi with external magnetic 
field h cx e z and parameter e, i.e., 



£ GL 

evaluated on the configuration 









2 


u,A' 




(v - iA'^j u 


+ 











curlA' — h c 



+ £ - 2 (i-N 2 ) 2 }, 



(u , A') = (u , h C x(e)e z A x/2) 



(5.33) 
(5.34) 



and (|5.30[) corresponds to the GL regime where the external magnetic field is between the first and the 
second critical fields. We can thus apply the lower bound for the GL functional proven in [2D], Theorem 
1.1 (note that in the definition of the GL functional given in [2D] there is overall factor 1/2), to get 



M > (1 - 0(1))^ log ' = (1 - O (l))0£ log > (1 _ o(1)) 



(5.35) 



since p TF (r, ; ) > u>\ log <5 1 ^ 1 inside % if Q > e -1 , and p TF (ri) > C, if SI < 



□ 



The proof of the lower bound to the GP energy E GP is now almost complete. Collecting the lower bounds 
inside all cells proven in the proposition above, we have 



£ GP \u] > 



m 2 1 log 7 1 



]Tp TF (rO(l- (l)). 



(5.36) 



fi£C 



The replacement of the Riemann sum by the integral can be done exactly as in (|5.18p : By the symmetry 
of the lattice cell and the L 1 — normalization of p TF , 



l-o(l) 

i 1 



(5.37) 



and we finally obtain 



Since 7 = min[£,e 2 f2] this gives the lower bounds in ()3.3p and 
l/(e 2 |loge|) entered in (|5.9p . 



E GP > E TF + ^f^-(l-o(l)). (5.38) 

Note that the condition il <C 



6 Vorticity of GP Minimizers 

In this Section we prove Theorem l3.3| which is a consequence of the energy asymptotics in (|3.3p together 
with a similar result in GL theory. Indeed we shall prove that one can associate to any GP minimizer 
a GL configuration satisfying certain energy bounds, which, by exploiting a result proven in [20] . yield 
the uniform distribution of vorticity. The proof closely follows the analysis performed in Section 5 in [20j 
and relies on Proposition 5.1 in this reference as a key ingredient. 

It is appropriate to point out that Theorem 13.31 is a statement about the uniform distribution of the 
local winding numbers of \& in the support of p TF but not about the nature of the singularties. The energy 
considerations behind this result are not sufficient to exclude the occurrence of singularities that are not 
pointlike, e.g., lines of zeros of ^/ GP . While we expect that v]/ GP contains only isolated vortices in the 
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parameter range (jl.ip . a proof of this has not been accomplished. The same is true for the corresponding 
question in GL theory (see, e.g., [2"TJI |2"T] ). 

Proof of Theorem EOl 

Exploiting the energy bounds proved before (see (|4.49j) . ()5.15|) . ()5.18j) and (I5.22p ). we obtain, with £ and 
t as in (pHl)| - (p3TT]) . 



o(l)^ 2 |log(e 2 ft)| 



TF 



(n)< 



TF 



?i£C 



(n) 



f^ 2 |log( £ 2 ^)l 



< cn, 



(6.1) 



so that, since the sum is performed over the lattice C C T, i.e., where p TF > cu\ log J| 1 , 



^ 2 |log(e 2 r!)| 



<g(e) ni 2 \\og(e 2 n)\J2p TF (n 



(6.2) 



?i£C 



for some g(e) — > 0, as e — * 0. 

Now we can distinguish, as in [20] . between good and bad cells, where the above inequality yields an 
upper (resp. lower) bound. The key point is that, if the definition of such cells is done in the appropriate 
way, the upper bound can be used to prove a uniform distribution of vorticity (inside good cells) and, 
at the same time, there are only few bad cells, i.e., their number is only a remainder with respect to the 
total number of cells. The final result would then be a simple consequence of the fact that cells cover 
supp (/Q TF ) in the limit e — > 0. 
We say that a cell Q l is a good cell, if 



(6.3) 



while inside bad cells the inequality is reversed. 

We can thus apply to any good cell Proposition 5.1 in [5D], which implies the existence of a finite family 
of disjoint discs B % e , i = 1, . . . , k, such that the sum of all the radii is bounded by f2 -1 / 2 and \u\ > 1/2 on 



dB\. Points 1, 2 and 3 in Proposition 13.31 then easily follows. In particular point 2 follows from a simple 
bound on the total number of cells, i.e., N <C which is a consequence of the conditions (|5.17p on £. 
Furthermore, setting d l equal to the winding number of u on dB\, which is also the winding number of 
j^GPj-i^GP Decause ^tf > q ms ^ e T, we have 



2tt2^<£ > fir(i-o(i)), 



27t2J4| <or(i + o(i)). 



(6.4) 



The second estimate above in particular implies that, by (|5.17|) . the measure on the left hand side of 
(|3.5j) is uniformly bounded in e, which guarantees its weak convergence. It remains only to show that it 
converges to the uniform measure on supp (p TF ) ■ To this purpose we first have to show that the number 
of bad cells included in any given open set S C supp (p TF ) (independent of e, i.e., such that \S\ > C > 0) 
is, for e — ► 0, much smaller than the total number of cells in S. In fact, since the area of S is positive and 
the diameter of the cells tends to zero for e — > 0, it suffices to shows that this is true for S = supp (p TF ). 
Denote by X the sets of indices i £ N such that Q l is a good (resp. bad) cell, then, by definition of bad 
cell and |Q) , 



N B nt 



|io g ( £ 2 tt)|V^)<£ 



^ 2 |log(e 2 fi)| 



< CNnl 2 \log(e 2 n)\g(e) 



(6.5) 
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where N B denotes the number of bad cells and N the total number of cells. As a consequence 

^V B < C^(e)N, (6.6) 

and, since g(e) = o(l), the number of bad cells is always much smaller than the total number of cells. 
The result can be easily extended to any set S C B\ by observing that the upper and lower bounds to 
the energy applies to any open subset of B\ . 

Theqrefore, for any given open subset S C supp (p TF ), good cells exhaust the whole of S as e — ► 0, i.e., 
N G £ 2 — > |iS|. Now, collecting all the disc families inside good cells and setting 

V= —^2 d i,e5(r-ri,s) , (6.7) 
where r^ e stands for the center of B l £ , one has, by the first estimate in (|6.4|) . 

H(S)>N G £ 2 > (l-o(l))|5|, (6.8) 
and similarly, by the second estimate in (16.41) . 

fi(S) < (1 + o(l))N G £ 2 < (1 + o(l))|5|, (6.9) 
which implies (|3.5| . since 5 is arbitrary. 

□ 



7 Conclusions 

Within the framework of two-dimensional GP theory we have evaluated exactly to subleading order the 
contributions of vorticity to the energy of a rapidly rotating Bose-Einstein condensate in a finite, flat 
trap. The results of the mathematical analysis lend support to the physical picture of a large number of 
vortices that are arranged in a triangular lattice at not too high rotational velocities but are eventually 
replaced by a 'giant vortex' with all the vorticity located outside the bulk of the density at sufficiently 
fast rotation. It would be desirable to substantiate this picture even further by generalizing Theorem 
13.31 for fl 3> £ _1 and by proving a lower bound to the energy matching the 'giant vortex' upper bound 
of [TB] to subleading order for f2 > e~ 2 \ loge| _1 . Further interesting open problems concern the nature 
of the singularities of the GP minimizer, in particular the exclusion of line singularities and the precise 
arrangement of the vortices. This would in particular require energy estimates beyond the subleading 
order considered here. 
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A The TF Energy and Density 

We collect here from [16] some formulas for the TF energy and density (in a slightly different notation). 
Defining 

u h = (A.l) 
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we have 



e 2 E TF = 



1 UJ Z TTUJ 

, II UJ < UJh, 

7T 8 768' ~ ' 



(A.2) 



if uj > uj\i, 



p TF (r) = < 



1 UJ 2 UJ 1 , . 



(A.3) 



, if uj > uj h , 



20F 8 

where [t]+ = t, if t > 0, and otherwise. The TF density p TF can be as well expressed as 

,2^.2" 



(A.4) 



where the chemical potential /i = i£ + e II/ 9 II2 i s nxe d by the normalization of p and it is 
explicitly given by 



2_lu 2 

_uS 
4 



1 - 



if uj < uj^ 



if uj > ujh- 



(A.5) 



Note that, if uj > uj^, a 'hole' centered at the origin occurs in the TF minimizer, i.e., p TF (r) = for 
all r < i?h with 

*«(»-sr <"> 

the radius of the hole. For uj > ujh and i?h < r < 1 we can also write the density as 



(A.7) 



Note also that the behaviour of E TF in the regimes f2 <C £ 1 (uj — > 0) and f2 3> e 1 (w — > oo) is 
respectively 



2 t^TF 



£ 2 £ 



0(w 2 ), fH£ 



2 cnTF 



or 



(l-O^- 1 )). 



(A. 
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